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ABSTRACT 

The emission from a gamma-ray burst (GRB) photosphere can give rise to a va- 
riety of spectral shapes. The spectrum can retain the shape of a Planck function or it 
can be broadened and have the shape of a Band function. This fact is best illustrated 
by studying GRB090902B: The main 7-ray spectral component is initially close to a 
Planck function, which can only be explained by emission from the jet photosphere. 
Later, the same component evolves into a broader Band function. This burst thus 
provides observational evidence that the photosphere can give rise to a non-thermal 
spectrum. We show that such a broadening is most naturally explained by subpho- 
tospheric dissipation in the jet. The broadening mainly depends on the strength and 
location of the dissipation, on the magnetic field strength, and on the relation be- 
tween the energy densities of thermal photons and of the electrons. We suggest that 
the evolution in spectral shape observed in GRB090902B is due to a decrease of the 
bulk Lorentz factor of the flow, leading to the main dissipation becoming subphoto- 
spheric. Such a change in the flow parameters can also explain the correlation observed 
between the peak energy of the spectrum and low-energy power law slope, a, a cor- 
relation commonly observed in GRBs. We conclude that photospheric emission could 
indeed be a ubiquitous feature during the prompt phase in GRBs and play a decisive 
role in creating the diverse spectral shapes and spectral evolutions that are observed. 

Key words: gamma-rays: bursts - radiation mechanism: thermal 



1 INTRODUCTION 

The original fireball model of gamma-ray bursts (GRBs) pre- 
dicts a strong photospheric component during the prompt 
phase (Goodman 1986, Paczynski 1986). The very high opti- 
cal depth to scattering expected near the base of the flow im- 
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plies that, regardless of the exact nature of the emission pro- 
cess, the resulting spectrum thermalises and is observed as a 
Planck spectrum. However, only a few GRBs have been iden- 
tified to be dominated by a Planck spectrum (Ryde 2004). 
Non-thermal spectra are more typically observed (Preece et 
al. 1998, Kaneko et al. 2006). Moreover, if the photosphere 
occurs far from where the acceleration of the flow ceases (the 
saturation radius), the thermal component is weakened by 
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adiabatic expansion. Most of the flow energy is then in the 
form of kinetic energy and only a thermal relic is left. 

It was therefore argued that the dominating emission 
mechanism should instead be optically-thin synchrotron- 
emission (Tavani 1996), emitted by relativistic electrons. 
These are accelerated following kinetic energy dissipation 
that takes place above the photosphere. Such dissipation 
could, for instance, result from internal shocks within the 
flow (e.g., Rees & Meszaros 1994). However, this paradigm 
has, in its turn, several severe problems. Foremost, in or- 
der to reproduce the observed sub-MeV spectral peak by 
synchrotron emission, a strong magnetic field, typically of 
the order B ~ 10^ — 10^ G is required. In such a strong 
magnetic field, the shocked electron population is expected 
to cool rapidly, producing a typical spectrum which is in 
stark contradiction to the observed spectral shape (Crider 
et al. 1997, Preece et al. 1998, Ghiselhni & Celotti (1999); 
see further discussion in e.g. Ryde et al. 2006). 

The challenges for optically-thin synchrotron emission 
led to the revival of the idea that the jet photosphere may 
play an important role, in one way or another, in the forma- 
tion of the spectrum (Eichler & Levinson 2000, Meszaros & 
Rees, 2000). Meszaros & Rees (2000) proposed that several 
spectral components exist in the gamma-ray band, includ- 
ing a Planck spectrum from the photosphere. Such a spec- 
trum can appear as a Band function when it is observed in 
a narrow energy band. Indeed, Ryde (2005), who studied 
time-resolved spectra from subpulses, showed that in many 
cases the GRB spectra in the BATSE energy range (~ 25 
- 1900 keV) are statistically indistinguishable between fits 
with a Band function and with a Planck function combined 
with a power- law (BB-|-pl). In many cases BB4-pl model is 
even preferred over the Band function. In addition, the ther- 
mal component was found, in these fits, to have a recurring 
behaviour during individual pulses: the temperature decay 
is well fitted by a characteristic broken power-law in time 
(see further Ryde & Pe'er, 2009). Even though the BATSE 
observations were limited by the narrow energy band, they 
thus gave an indication that a photospheric emission does, 
in fact, exist in many bursts, and that GRBs in general have 
several spectral components in the gamma-ray band. 

Moreover, Ryde & Pe'er (2009) argued that GRB spec- 
tra should typically be more complicated than a single Band 
function, when observed over a broader energy range. This 
is because approximately 10% of all CGRO BATSE bursts 
have, during their entire duration, time-resolved spectra 
with a Band function high-energy spectral-slope /3 > —2 
(Kaneko et al. 2006). The peak in the energy flux iyF^) of 
these bursts must therefore be above the BATSE energy 
range, at an energy higher than the determined spectral 
break. Ryde & Pe'er (2009) also studied the few bursts for 
which there are simultaneous and time-resolved data avail- 
able from both BATSE and EGRET- TASC. They found 
that several breaks indeed exist in the spectrum (see also 
Barat et al. 1998). While the overall power peak, in the stud- 
ied bursts, lay in the EGRET range, Ryde & Pe'er (2009) 
interpreted the break in the BATSE range as a subdominant 
thermal peak. Recently, similar conclusions were drawn by 
Guiriec et al. (2011) and Zhang et al. (2010). 

There are thus strong arguments, both theoretical and 
observational, that the photospheric emission plays an im- 
portant role in the spectral formation during the prompt 



phase in GRBs. Indeed, recently the Fermi Gamma-ray 
Space Telescope (energy range 8 keV — > 300 GeV) has ob- 
served bursts in which a photospheric component is present 
and several spectral components are required (Ryde et al. 
2010, Guiriec et al. 2011). However, many bursts observed by 
the Fermi instruments do not show such obvious distinction 
between a Planck function and other spectral components. 
The spectral peak is more typically described by a more 
broadly peaked Band component. 

This has, in its turn, led to increased interest in the 
suggestion put forward by Rees & Meszaros (2005), who 
argued that strong dissipation should naturally occur below 
and close to the photosphere (see further ij3|. This results 
from oblique shock waves, that are formed at the edges of the 
jet as it propagates through the star. The dissipated energy 
results in reprocessing of the original Planck spectrum, due 
to, e.g., Comptonization by energetic electrons (Pe'er et al. 
2005, 2006; see also Giannios 2006). This alters the Planck 
spectrum into the observed spectrum. Emission from the 
photosphere could therefore have a Band-like character and 
its shape should depend on the details of the dissipation of 
the kinetic energy of the flow. Similar scenarios have recently 
been discussed by several authors, e.g., Beloborodov et al. 
(2010), loka (2010), Lazzati & Begelman (2010), Toma et al. 
(2010), and Bromberg et al. (2011); see further discussions 
in Ramirez-Ruiz (2005) and Ruffini et al. (2005). 

In this paper, we study the particular spectral evolution 
of GRB090902B. In this burst, the main spectral compo- 
nent, stemming from the photosphere, exhibits a change in 
spectral characteristics half-way through the prompt phase. 
At early times, the thermal component resembles a Planck 
function, while at late times this component broadens sig- 
niflcantly. Based on this study, we further discuss the condi- 
tions under which photospheric emission is broadened. We 
argue that this mechanism can be applicable to more typi- 
cal spectral evolutions, in which the spectra more gradually 
evolve from being hard to becoming softer. In particular, we 
argue that it may provide a natural explanation for the ob- 
served variety of spectral shapes in GRBs (in particular the 
width of the spectral peak). In §2 the spectral behavior of 
GRB090902B is presented and in §3 the effects of subphoto- 
spheric shocks on GRB spectra are discussed. In §4 we use 
subphotospheric shocks to explain the spectral evolution in 
GRB090902B. We discuss our results in §5, and we conclude 
in §6. 



2 SPECTRAL BEHAVIOR OF GRB090902B 

The bright and long burst GRB090902B was detected by the 
Fermi Gamma-ray Space Telescope by its two instruments 

- the Large Area Telescope (LAT; energy range 100 MeV 

- > 300 GeV) and the Gamma-ray Burst Monitor (GBM; 
8 keV - 40 MeV). The burst lies at a redshift of z = 1.822 
(Cucchiara et al. 2009). It is one of the strongest bursts 
detected by Fermi and the emission at energies larger than 
8 keV lasted for approximately 25 s. The most energetic 
photon with an energy 

of 33.4j:^-7 GeV was detected at 82 s 
after the trigger by the LAT. The light curve of the prompt 
phase is shown in Figure [l] 

During the prompt phase of approximately 25 s two 
distinct, separate components are observed throughout the 
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Figure 1. Light curve of GRB090902B, observed by Nal 1 in the 
energy range 8.5-904 keV. The presented time-binning is used in 
the analysis and the two epochs discussed in the paper are divided 
by the vertical line. The thin dashed line is the background level. 



duration: a peaked MeV component (modeled with a Band 
function) and a power-law component. The power-law com- 
ponent is clearly detected at energies both below and above 
the MeV peak (observations are made in the range of 8 
keV - ~ 30 GeV). Moreover, while the MeV peak undergoes 
substantial spectral evolution the power-law component re- 
mains relatively steady with photon-index of approximately 
-1.9 (Abdo et al. 2009). 

At early times (first 12.5 s; epoch 1) the pronounced 
MeV peak is so steep and narrow, that it must be attributed 
to emission from the photosphere (e.g. Fig. 1 in Ryde et 
al. 2010, see also {2.2 I. However, during the second half of 



the prompt phase, which lasted for another 12 s (epoch 2), 
the MeV peak differs significantly from a Planck function, 
resembling a typical Band spectruirj^ 

We argue here that, although the spectral shape of the 
MeV peak varies, it is most likely that it has the same origin 
throughout the burst. The reasons are the following, (i) The 
spectra are clearly separated into two spectral components 
throughout the burst duration, namely an MeV peak and an 
independent power-law component (Fig. |3|. (ii) Although 
the MeV peak broadens, the power-law component remains 
relatively steady (see Abdo et al. 2009). (iii) The spectra 
of the MeV bump during epoch 2 are still inconsistent with 
the expected non-thermal spectrum; fast cooling electrons 
yields a = —1.5 (see [2.2 below), (iv) To get a synchrotron 



peak-energy to lie in a similar energy range as the thermal 
peak during epoch 1 requires an unreasonable coincidence, 
(v) On the other hand, as we show in ^ below, broadening 
of the thermal peak by sub-photospheric dissipation easily 
reproduces the observed spectrum. Therefore, we argue that 
this burst provides an observational evidence that the emis- 
sion from a GRB photosphere does not necessarily need to 
be a narrow, Planck-like spectrum but can be significantly 



broadened. This burst is the best example available to study 
the details of the photosphere and its emission. 

In general, the spectral evolution in GRBs is more grad- 
ual than observed in GRB090902B. The steepest sub-peak 
slopes (largest values of a) are typically found at the very 
beginning of the prompt phase, and only during a small frac- 
tion of the burst duration. The spectra thereafter rapidly 
soften (Crider et al. 1997, Ghirlanda et al. 2003). The par- 
ticular property of the spectral evolution in GRB090902B is 
the substantial fraction of the burst duration during which 
the emission spectrum is Planck-like. This allows the estab- 
lishment of its photospheric origin. 

Here, we further study the time-resolved spectra by fol- 
lowing the analysis performed in Ryde et al. (2010). We use 
the same detectors (Nal 0, 1 and BGO 0,1 and LAT front 
and back) and the time binning was determined by requiring 
a signal-to-noise ratio in the strongest illuminated detector 
to be at least SNR = 45. The data were fit using RMFFlQ 
version 3.0 using the Castor C-statistic to determine the 
goodness-of-fit. The light curve with the time binning used 
is shown in Figure [l] Here the burst is divided into two 
epochs: one comprising data from within 12.5 s of the trig- 
ger, referred to as epoch 1 (analyzed in Ryde et al. 2010), 
and another comprising data after 12.5 s until the end of the 
burst, referred to as epoch 2. 



2.1 Band function fits 

A model consisting of a Band function (Band et al. 1993) 
(for the dominating MeV component) and a power-law com- 
ponent fits the time-resolved spectra well during epoch 1 
(Ryde et al. 2010). This is also the case for our fits. How- 
ever, inspection of the C-stat maps of the parameters reveals 
that the error ranges are not well constrained in several bins. 
The reason for this was determined to be the low amplitude 
of the power-law in these bins. In these cases, we froze the 
power-law amplitude to the value found in the fit, and then 
performed a new fit to determine values and uncertainty 
intervals for the remaining parameters. For all bins, the val- 
ues of the parameters found before freezing the power-law 
amplitude were consistent with the values found after the 
amplitude was frozen. 

Table 1 shows the results of our fits to the time-resolved 
spectra indicating a (photon-index of the sub-peak power- 
law), P (photon-index of the super-peak power-law), Ep 
(peak energy), and the photon- index of the power-law com- 
ponent. Figure [2] shows the evolution of the shape of the 
MeV component. The upper panels show the parameters a 
and P over the full duration of the burst (0-25 s). The in- 
dex of a = —2/3 expected from optically-thin synchrotron 
emission is indicated (Rybicki & Lightman 1979) as well as 
a — 0, the spectral slope expected from, e.g. jitter radiation 
and from extreme cases of inverse Compton emission from 
a delta function distribution of electrons (Jones 1968). The 
errors on the data points represent asymmetric, one-sigma 
uncertainties on the parameter values found from the fitting. 

In Figure [2j it is clearly seen that the spectral shape 



^ Note that a Planck spectrum can in principle be approximated 
by a "Band" function, but with very steep a = 1 and f} — oo. 
We use the term "Band" to describe spectrum that is not as steep. 



2 R. S. Mallozzi, R. D. Preece, & M. S. Briggs, "RMFIT, A 
Lightcurve and Spectral Analysis Tool," ©Robert D. Preece, Uni- 
versity of Alabama in Huntsville. 
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Figure 2. Evolution of ttie MeV component in GRB090902B. Upper left panel. Evolution of low-energy photon-index a. Two horisontal 
lines are shown, which correspond to a = 0, the most extreme value expected for inverse Compton models and a = —2/3, expected for 
optically-thin synchrotron emission for a slow cooling electron population. The dashed line indicates 12.5 s. Upper right panel. Evolution 
of the high-energy power-law index /3. Lower left panel. Correlation between a and /3. Note that some points have only one-sided error 
bars indicating that they are unconstrained in the other direction. Lower right panel. Peak-aligned Band functions corresponding to the 
fits at two different times, illustrating the spectral broadening. These two times include the narrowest and the broadest Band spectra. 



exhibits a change in character at approximately 12.5 s after 
the trigger, as noted already by Ryde et al. (2010). From 
being very peaked, with a ~ 0.3 and /3 ~ —3.5 and with a 
spectral width of the pealj^w = 6, the spectrum broadens 
significantly. During the second epoch, the typical value of 
a ~ —0.6 and P ~ —2.5. As a result, the spectral width 
has increased between the first and second epochs by more 
than a factor of two, to typically w =10 - 20. In the lower 
right-hand panel in Figure [2] the best-fit Band function of 
the MeV peak is plotted for two instants, at 6.5 s and at 22 
s after the GBM trigger. These spectra have been aligned 
to each other's Ep values in order to highlight the spectral 
broadening: while the typical spectrum at epoch 1 is close to 



^ We here define the spectral width, w, as the ratio i^high/^'lowi 
where energy fluxes Fe{Ei^^) and -^^{-Bhigh) a-re equal to half 
the peak flux (-Fpoak/2) below and beyond the peak, respectively. 



a Planck function, the spectrum at epoch 2 has a shape that 
is typical for a GRB spectrum, that is, a Band function. 

Even though the change in spectral character between 
the epochs is clear from the above discussion, there are still 
important similarities between the epochs. Most notably, 
Abdo et al. (2009) clearly showed that the power-law compo- 
nent remains relatively steady during the spectral evolution 
of the prompt phase. Moreover, during the time period 15 
- 17 s, the MeV component becomes relatively hard again. 
The similarities between these spectra and the epoch 1 spec- 
tra are shown in Figure |3] Here a spectrum from epoch 1 
(8.1- 8.5 s) and a spectrum from epoch 2 (15.9 - 16.4 s) are 
shown. The broadening of the MeV component, compared 
to the epoch 1 spectra, is still apparent, even though the 
broadening is not as large as for the other epoch 2 spectra. 
These observations are thus strong indications of that the 
emission during the two epochs are of similar origin, i.e. a 
photospheric and an optically-thin component. 
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Figure 3. Time resolved uF^, spectrum for two time intervals t = 8.1 — 8.5 s (epoch 1) and t = 15.9 — 16.4 s (epoch 2). The Band + 
power-law model is fitted to the data over the GBM + LAT energy ranges. The symbols refer to the different instruments on the Fermi 
Gamma-ray Space Telescope. While there still are similarities between the spectra, the broadening of the MeV component is apparent 
(Compare to Fig. 1 in Ryde et al. 2010). 
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Figure 4. The evolution of the peak energy of the Band function 
fits, Ep, as a function of time since the GBM trigger. 



Figure [1] shows the Ep evolution as a function of time. 
The averaged value of Ep is lower during epoch 2 compared 
to epoch 1. 

In spite of the variations in the spectrum of the MeV 
component, it dominates the spectral energy flux throughout 
the burst duration. The ratio of the energy flux in the MeV 
peak relative to the total Fermi 7-ray flux (MeV peak -|- 
power law component) is in the range 80 and 95 % during the 
entire burst duration. This mainly results from a decrease 
of the amplitude of the PL component towards the end. 



2.2 Significance of the hard o-values. 

During epoch 1 the averaged value of a is < a >= 0.11 and 
for several of the time-bins the value of a is even steeper, 
with a ^ 0.2. These very hard spectra are particular chal- 
lenging for non-thermal models. Synchrotron and inverse 
Compton emission in the fast cooling regime are expected 
to produce a spectral slope of a = —1.5. Here, we consider 



the most extreme scenario with a = (Jones 1968, Epstein 
& Petrosian 1973). We therefore want to estimate the sig- 
nificance of rejecting a = for these spectra. This allows us 
to determine the significance of the conclusion that neither 
synchrotron nor inverse Compton scattering processes can 
explain the spectra. 

This is done by simulating spectra with RMFIT v3.0 
using the set of observed spectral parameter values of the 
Band function fit. We use the Fermi detector responses for 
GRB090902B and take into account the Poissonian nature of 
the observed counts and realistic background emission. The 
simulated count distributions are then fitted in the same 
way as the real data. The parameter values we find can then 
be studied. 

We illustrate the procedure on the 7th time-bin (6.3- 
7.0 s) which has the hardest observed value of a = 0.3 ± 0.1 
(Ryde et al. 2010). We start with the null hypothesis that 
the spectrum observed in this time bin has an actual value 
Q = 0. We therefore freeze a at this value, and, based on 
the best fit parameters we find from the data, we perform 
100,000 simulations. This large number of simulations allows 
us to make a proper estimation of the significance level. The 
simulated spectra are then fitted with the Band function 
with all parameters free to vary. The distribution of the val- 
ues of a that we find are shown in Figure [5] The distribution 
is approximately Gaussian and is slightly skewed. It has a 
mean value of a — 0.01 and a standard deviation of 0.068. 
The inset in Figure [5] is a magnification of the distribution 
around a = 0.3. This figure shows that 8 out of the 100,000 
simulated spectra have fits with values of a greater than 0.3. 
This gives the probability that a spectrum with an actual 
value a = should be observed with a > 0.3 by chance to 
be 8 X 10~^. Therefore the null hypothesis can be rejected 
on a very high confidence level for this bin. Most of the 
measured values of a during epoch 1 are, though, not signif- 
icantly inconsistent with a = (see Fig.[2|. In any case, the 
epoch 1 spectra are challenging for purely non-thermal emis- 
sion models, since a = is only expected under somewhat 
extreme conditions, see further |5.1[ The expected value is 
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Table 1. Results of spectral fits to the data of GRB090902B during the time 12.54-24.58 s after the burst trigger (epoch 2). Uncertainties 
marked unc indicate that the parameter is unconstrained. Results for epoch 1 are presented in Ryde et al. (2010). 

rather a — —1.5, which is produced by a population of fast 
cooling electrons (Ghisellini & Celotti, 1999). 

The averaged value of a during epoch 2 is < a > = 
—0.65. This value is still significantly inconsistent with a = 
— 1.5. For instance, the 42nd time-bin (19.84-20.22 s) has one 
of the softest values of a with a = — 0.86±0.04. Even though, 
when we perform a fit with a frozen at -1.5 to this time-bin, 
we find an increase of the C-stat value from 488 to 1167 for 
597 degrees of freedom. The hypothesis of having a = —1.5 
can therefore be rejected on a significance level of less than 
10"^". We also point at the fact that several of the fitted 
spectra in epoch 2 have a ~ —0.2; the 32nd time-bin (15.87- 
16.38 s) has the hardest spectrum with a = —0.2 ± 0.1, 
approaching the values of a from epoch 1, see Figure 2. 



3 SUBPHOTOSPHERIC HEATING AND ITS 
EFFECT ON THE EMISSION SPECTRUM 

As we saw above, GRB090902B is a particularly interesting 
burst since the main spectral component initially is close 
to a Planck function but later evolves into a broader Band 
function. Indeed, as we will see in this section, the emitted 
spectrum that is expected from a GRB photosphere depends 
on the existence of significant dissipation close to the pho- 
tosphere. 




Figure 5. Histogram distribution of a values found from 100,000 
simulated spectra. The inlay is a magnification of an interval 
around a = 0.3. See the text for details. 



3.1 Photospheres and dissipation radius 

The photosphere related to the electrons associated with 
the baryons in the outfiow (the baryonic photosphere) is 
expected to be at (e.g. Meszaros et al. 2002) 



rph ~ 4.8 X 10' 



L 



1053erK/s (63o) 



(1) 
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where the typical value for GRB090902B is used: L = 10^^ 
erg/s (Abdo et al. 2009). We also use the time- averaged value 
for r = 600 F^/'' ~ 630 (Ryde et al. 201of] Here Y de- 
notes the ratio between the total fireball energy and the 
energy emitted in 7-rays and can be estimated from after- 
glow measurements. Cenko et al. (2010) estimated the value 
to r ~ 1.2 in the case of GRB090902B. 

Pair formation can be significant and modify the spec- 
trum (e.g. Eichler & Levinson 2000, Pe'er & Waxman 2004). 
Subphotospheric dissipation can lead to copious production 
of pairs, resulting in a second, pair photosphere above the 
baryonic photosphere (Meszaros et al. 2002). During the dis- 
sipation process, electrons are expected to be accelerated to 
high energies, thereby emitting energetic photons at ener- 
gies above the threshold for pair production, me(? ■ These 
photons then produce pairs by annihilating with the lower 
energy photons. The created pairs have modest Lorentz fac- 
tors, 7e ~ few (in the comoving frame). For sub-photospheric 
dissipation these pairs are expected to be more numerous 
than the baryon-related electrons, and a pair photosphere 
is expected to be established at radius (Pe'er & Waxman, 
2004) 



V630y 



(3) 



r± ~ 8.0 X 10" 



L 

lO^-'crg/s 



t± a 



(63o) 



(2) 



where e± is the fraction of the total fireball energy that is dis- 
sipated into photons with comoving energy larger than m^c? 
and is available for pair formation. We define ro = ar^ as 
the effective radius at which the outflow starts to accelerate 
(in the absence of dissipation): a ^1 and Tg = 2GM/c? — 
3 X 1O'^(M/1OM0) cm is the Schwarzschild radius of the 
central black hole of mass M — 10 Mq. The pair photo- 
sphere is above the baryonic photosphere if e± > rric/mp — 
5.45 X 10"'' (Rees & Meszaros 2005). 

The existence of a dissipation process is in fact required 
by the data. Since we see non-thermal spectra in GRBs this 
implies that there is a mechanism that dissipates some frac- 
tion of the jet kinetic energy. The exact nature of this dis- 
sipation is debatable. Several models exist in the literature. 
The leading ones are: (i) internal shocks in which shells with 
varying Lorentz factors interact with each other (Rees & 
Meszaros 1994); (ii) oblique shocks within the funnel in the 
star (e.g. Morsony et al. 2007); (iii) coUisional dissipation in 
the flow (Beloborodov 2003); or (iv) in the case of Poynting- 
flux dominated flows, the dissipation could result from mag- 
netic reconnection (e.g. Thompson 1994, Giannios & Spruit 
2005). The shocked region subsequently cools by emitting 
photons through, for instance, synchrotron emission and/or 
inverse Compton emission. 

In the internal shocks scenario, dissipative heating is 
typically assumed to occur well above the photosphere: with 
variations in Lorentz factor in the flow of the size AP ~ P 
the internal shocks occur at r^h ~ 2ror^ ~ 2.1 x lO" cm. 
However, Rees & Meszaros (2005) pointed out that due to 
jet edge effects oblique shocks might form below and close 
to the photosphere: 



Later we also use the notation r2.8 = (r/630) and L53 = 
(L/1053erg/s): Note that Pe'er et al. (2010) used F ~ 1000 which 
was estimated for a single time bin used in their analysis 



3 X 10-2 

where 9j is the nozzle opening half angle. 

It thus follows that internal shocks may naturally be 
expected in proximity of the photosphere. While the details 
of these process are still highly uncertain, this problem is ex- 
tensively being studied numerically. Indeed, numerical simu- 
lations of a jet penetrating though the core of the progenitor 
show that such shocks do indeed occur (Lazzati et al. 2009, 
Mizuta et al. 2010). We note that similar conclusions are 
drawn for other dissipation processes as well (Giannios & 
Spruit 2005, Giannios 2006, 2007, Tchekhovskoy et al. 2008, 
Beloborodov 2010). 

3.2 Broadening of the Planck spectrum: 
analytical arguments 

The energy that is dissipated in subphotospheric shocks 
partly thermalises again to an extent that depends on the 
conditions at the dissipation site, particularly on the optical 
depth. Detailed calculations of the thermalisation process of 
such shocks by Pe'er et al. (2006), and in S3.3 below show 



that a large variety of spectral shapes can be achieved. Sim- 
ilar results are found in magnetic reconnection models (e.g., 
Giannios 2006). In particular, Pe'er et al. (2006) showed 
that the Planck function that is injected into the dissipa- 
tion region is modified to a varying extent depending on 
the dissipation fractions and the optical depth. The Planck 
spectrum therefore loses its original shape and the outgoing 
photospheric emission has a non-thermal shape. The result- 
ing spectrum can have a rather complex spectral shape. As 
we show here, under plausible conditions in many cases it 
can be described as a smoothly broken power law. 

While detailed numerical results are presented below, 
we give here some basic analytical arguments to describe 
the conditions under which significant modification of the 
spectrum can take place. 

Assume that the dissipation process, regardless of its ex- 
act nature, produces a population of energetic electrons with 
characteristic Lorentz factor 7^ 3> 1. These electrons cool by 
Compton scattering the thermal (photospheric) photons, on 
a time scale given by tcooi — mec/ {4:/3)yeaTUph- Here, Uph 
is the energy density in the thermal photon component, and 
(jt is Thomson cross section. This loss time can be com- 
pared to the dynamical time of the problem, tdyn = r/Tc to 
obtain (see Pe'er et al. 2005) 

tloss _ 3 TTteC^r _ 3 Uel 
tdyn A'yeOTUphr 1: ^lUphT^e' 

Here, r is the dissipation radius, r-^e ~ {r /T)neOT is the 
optical depth to photon scattering by the electrons and 
Uel ~ 'yeneTrieC^ is the energy density in the electron compo- 
nent. Equation Q thus implies that, regardless of the exact 
nature of the dissipation process, the electrons cool rapidly, 
on a time scale much shorter than the dynamical time. This 
is provided that (i) the dissipation occurs below or not too 
high above the photosphere, and (ii) that the energy that is 
being dissipated is not much larger than the energy density 
in the photosphere (see Pe'er et al. 2005). The rapid cooling 
is due to Comptonisation of the thermal component; obvi- 
ously, synchrotron emission further contributes to the rapid 
cooling of the electrons. 
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This rapid cooling implies that the electrons lose most 
of the energy imparted to them by the dissipation process. 
This energy is used to up-scatter the photospheric photons, 
as well as to emit synchrotron photons. Part of this energy 
is converted into pairs, by upscattered photons that are en- 
ergetic enough. 

The distribution of the rapidly cooled electrons reaches 
a quasi-steady state: when the electrons are cold enough, 
inverse Compton scattering becomes inefficient, while other 
processes, such as direct Compton heating or synchrotron 
self-absorption heats the (cold) electrons. The electron 
distribution can therefore be approximated as a (quasi-) 
Maxwellian distribution, with characteristic temperature Te- 
As long as the dissipative process that heats the electrons 
(or introduces a population of energetic electrons into the 
plasma) exists, the steady temperature of the electrons is 
inevitably higher than the temperature of the photospheric 
(thermal) photons: Te ^T^. 

The plasma is therefore characterized by two tempera- 
tures, T-y and Te- Due to the rapid cooling, during most of 
the dynamical time the scatterings take place between the 
thermal photons and the cold electrons. Since ^Tj, the 
thermal photons gain energy, resulting in modification of the 
Wien part of the Planck spectrum, to produce a smoother 
cutoff at high energies (above the thermal peak) . The exact 
shape of the spectrum at these energies (which corresponds 
to in Band fits), depends on the optical depth, and the 
ratio of the energy densities in the electrons and photon com- 
ponents. A significant shift can be obtained if the Compton 
(or equivalently Kompaneet's) j/-parameter is of the order 
of a few (which translates to optical depth of a few), and 
energy densities in the electron and photon components are 
roughly comparable. 

Interestingly, somewhat similar conditions are required 
in order to obtain a significant modification to the Rayleigh- 
Jeans part of the Planck spectrum. This part is modified if 
two conditions are met: first, significant number of photons 
at energies below the thermal peak must be introduced into 
the plasma. The most natural way to obtain a large density 
of cold photons is via synchrotron radiation. This emission 
is expected at low energies, resulting from emission from the 
cold electrons. In order to obtain a significant flux, a strong 
magnetic field is thus needed. The second condition is that 
up-scattering of these photons leads to energies comparable 
to the original thermal photons. The condition here is again 
y ^ 1, which, due to the low value of Te is translated into 

T^e ~ few. 

We further note that if the optical depth T-,e — > <x, 
the spectrum approaches either a Planck or Wien spectrum, 
as the energy given to the electrons is distributed among 
the electrons and the photons. If the energy given to the 
electrons by the dissipation process is much larger than the 
energy in the thermal component, u^i ^ Uph, then the non- 
thermal part of the spectrum is significantly more luminous 
than the thermal part. In such a scenario, the thermal part 
may not be detectable (for very high optical depth, the re- 
sulting Planck or Wien spectrum will have a temperature 
that is different than the original temperature of the photo- 
sphere). 

We thus conclude that broadening of the Planck spec- 
trum naturally occurs if the following conditions are met: (i) 
Dissipation processes take place below the photosphere, at 



optical depth of r^e ~ few ; (ii) the energy given to the elec- 
trons is comparable to the energy in the thermal photons 
component; and (iii) a strong magnetic field exists, of the 
order us/iith ~ tens %. See further discussions in Pe'er et 
al. (2005, 2006), Giannios (2006), and Beloborodov (2010). 

3.3 Broadening of the Planck spectrum: Detailed 
numerical simulations 

The arguments given in |3.2| above provide a guideline to 
possible conditions that can lead to broadening of the ther- 
mal spectrum. However, quantitative results can only be 
obtained numerically. This is because of the non-linearity 
of the problem. First, a large number of pairs can in princi- 
ple be produced. Thus, a rapid electromagnetic cascade may 
be presented. Second, as most of the scatterings occur with 
cold electrons, the cross section is Klein-Nishina suppressed, 
and hence simple analytical approximations to the resulting 
spectra are absent. Finally, photons and electrons can par- 
ticipate in a large number of processes, such as synchrotron, 
synchrotron self-absorption and Compton scattering, which 
can have similar importance. 

In order to obtain numerical results, we use the code 
developed by Pe'er & Waxman (2005) for the study of GRB 
prompt emission. This code was further modified to the 
study of photospheric emission by Pe'er et al. (2005, 2006). 
The numerical code solves self-consistently the kinetic equa- 
tions that describe a large number of physical processes (syn- 
chrotron emission, synchrotron self-absorption, direct and 
inverse Compton scattering, pair production and annihila- 
tion and the development of an electromagnetic cascade) 
that can take place following the injection of energetic par- 
ticles close to the photosphere. Its unique integrator enables 
it to follow the evolution of the particle distribution and 
spectra over many orders of magnitude in time and hence 
in energy scales. Thus, it is ideal for studying processes that 
can take place in regions of high optical depth, in which 
the characteristic time scale for interactions is much shorter 
than the dynamical time. A full description of the code ap- 
pears in Pe'er & Waxman (2005). 

In the scenario considered in our calculations, the exact 
nature of the dissipation process is not specified. We assume 
that the outfiow is characterised by steady luminosity Lq, 
and constant Lorentz factor, F. The dissipation is assumed 
to take place at radius r^, and dissipate some fraction ed 
of the kinetic energy. A fraction ee of this energy is used 
to accelerate electrons while a fraction e_g is channeled into 
magnetic fields. 

The thermal component is considered as a constant 
source of thermal photons, which irradiate the interaction 
region during the whole calculation, and interact with the 
accelerated electrons. After the energy has been dissipated 
the evolution of the electron/positron populations and the 
photon population is followed during one dynamical time, 
tdyn ~ ri/Fc, subject to synchrotron radiation, Compton 
scattering and pair production and annihilation, as well as 
the constant influx of thermal photons. At the end of the 
dynamical time all the radiation is assumed to be released, 
and the emerging spectrum is given by the photon distribu- 
tion at that time. Since r ^ 1 at the dissipation site, one 
should take into account adiabatic expansion until the pho- 
tons are released at r ~ 1. However, since the optical depth 
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is not more than a few, adiabatic expansion only marginally 
affects the spectrum and can be neglected for our purposes 
(see further Pe'er & Waxman 2004). 

The code allows us to quantitatively confirm the analyt- 
ical results discussed in j ]3.2| For instance, while Compton- 
isation naturally leads to a harder high-energy power-law, 
we find that the optical depth has to be relatively high for 
the number of scatterings to be large enough to affect the 
high-energy slope. On the other hand, even larger values 
of T leads instead to a steepening of P due to thermalisa- 
tion. Moreover, the slopes of the low-energy power-law in- 
dex, a, depends most strongly (for a individual dissipation 
scenario) on the strength of the magnetic field generated, i.e. 
on e_B, giving rise to synchrotron emission. A larger es leads 
to a softer value of a. However, the effect of Comptonisa- 
tion and optical depth counteract this softening. Therefore, 
the broadest photospheric spectra are obtained for a strong 
dissipation occurring at moderate optical depths, typically 
r ~ 10. A more detailed accounting of the effects of subpho- 
tospheric heating on the photospheric spectrum, using this 
code, is given in Nymark et al. (2011). 

If shocks occur above the baryonic photosphere (r << 
1) we find that the original Planck function will only be 
marginally affected by the shocked electrons due to the low 
number of scatterings that will occur. 

Note that in these calculations we consider only one dis- 
sipation episode. In reality several dissipations are expected, 
making the emerging spectrum a superposition of the spec- 
tra from several dissipation episodes (see e.g. Giannios 2008, 
Pe'er et al. 2011). However, for the purposes of this study, 
the spectrum resulting from one dissipation episode is suffi- 
cient to get an indication of the effects that dissipation can 
have. 



APPLICATION TO THE SPECTRAL 
EVOLUTION IN GRB 090902B 



4.1 Epoch 1 



In the data analysis above ([2.1 \ we found that the averaged 
value of the low-energy slope is harder than a — during 
epoch 1 in GRB090902B. Apart from being a constraint on 
the radiation process, this observation sets constraints on 
any dissipation process that can have had an infiuence on 
the spectral shape, such as constraints on the dissipation 
fractions of the kinetic energy and dissipation radius. 

It is worth to note here, though, that when observing a 
GRB photosphere we actually expect to observe a superpo- 
sition of spectra with different temperatures as measured in 
the observer frame, due to geometrical effects (Pe'er 2008). 
The result is a slightly broadened spectrum. Therefore the 
sharpest spectrum is not a Planck function but rather a 
multi-color blackbody (Pe'er & Ryde 2011, see further Lund- 
man et al. 2011). This effect was proposed to explain the 
deviations from a Planck function that is observed during 
epoch 1 (Ryde et al. 2010, see also Larsson et al. 2010). 

In any case, the simulations of subphotospheric shocks 
and their effect on the emitted photospheric photon ener- 
gies, described in 33.3[ show that dissipation occurring at a 
low optical-depth only marginally affects the emitted pho- 
tospheric emission. Moreover, in the case of dissipation at 



moderate optical depths, say r = 10, a small deviation from 
a Planck function can only be achieved if the dissipation is 
relatively weak. Indeed, such dissipation can easily repro- 
duce the observed spectral shapes during epoch 1. As an 
example we show in the left-hand panel of Figure [6] a spec- 
trum found from simulating a shock dissipation in a flow 
at r = 10, and with the parameters = 0.1, — 0.1, 
e_B ~ 0.1. The peak in this spectrum is slightly broader than 
a Planck function (indicated in the figure with a red dashed 
line) . This shape of the spectrum indeed resembled the Me V 
peak observed in GRB090902B, e.g. see the red spectrum in 
the lower right panel of Figure [2j 

Since the observed spectra do not differ much from a 
Planck spectrum during this epoch, one can thus argue that 
there cannot have been significant subphotospheric dissipa- 
tion during this period. Only weak dissipation or no dis- 
sipation at all should have occurred. The thermal peak is 
therefore mainly the original non-processed thermal emis- 
sion, which was formed at (or close to) the base of the flow, 
at very high optical depths. 

The thermal emission that is observed can either be at- 
tributed to a baryonic photosphere or to a pair photosphere. 
Depending on which photosphere is assumed, a different con- 
clusion will be drawn. 



4.1.1 Baryonic photosphere 

The photospheric radius deduced from the observations to 
Tph ~ 10^^ cm can be assumed to be associated to the opac- 
ity of the baryonic electrons. This is what is assumed in Ryde 
et al. (2010) and Pe'er et al. (2010). Furthermore, since the 
photosphere can be assumed to be non-dissipative one can 
use the standard theory discussed in Pe'er et al (2007) and 
derive that ro = 10^ F"^''^cm ~ 7.6 x lO^cm (using the 
time average values for GRB090902B). The value of ro im- 
plies an 5 ~ 3.3 x 10^ Y'^^^M-^ ~ 2.5 x lO^Aff \ with the 
notation Mi = (M/10Mq). According to equation Q the 
large value of a, and thereby ro, implies that the dominant 
fraction of internal shocks occur well outside of the photo- 
sphere, assuming typical size of the opening angle 9j. This is 
consistent with the assumption of a non-dissipative outflow. 



4.1.2 Pair photosphere 

The value of the photospheric radius, ~ 10^^ cm, that is 
deduced from the observations can instead be attributed to 
the the pair photosphere. Thus setting r± ~ 10^^ cm then 
equation ^ yieldfj^ that e± a"^ = 1.25 x lO'^L^g^Fi g. 

Furthermore, equations ([T]) and ||2| imply that the ratio 
of the photospheres is r±/rph ~ 1.8 x 10'^ e± ci~^ in the 
coasting phase. Thus r±/rph = 2.1 L^^T 2. s and thereby the 



baryon-photospheric radius is rph 
Furthermore, ro = — 1 
saturation radius 1 



5.0 X W^'LsiT^ s cm. 
2.4 X We^MiLssT , ocm and the 



-- Fro = 1.5 X lO'^eiMiLssPa.gCm. 
Pairs are either created at the base of the flow (Good- 
man 1986) or in the energy dissipation at the shock. The 
requirement that rsh ~ r± (pairs are created in the shock), 
with rsh given by equation |3|, then corresponds to a ^ 1.6. 



Using r ~ 1000, which is inferred from opacity arguments 
(Abdo et al. 2009), rph ~ 2 X lO^^ cm and e± a'^ ~ 5 X 10"^ 
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From the estimate of e± a above this is equivalent of 
e± < 1.58 X 10~^. This value implies that in this scenario 
ro ~ 3 X lO^cm and ~ 2 x 10*^ cm for GRB090902B. 
This implies that the characteristic initial temperature of 
the fireball is of the order To ~ 4 MeV (e.g. Rees & 
Meszaros 2005). Therefore, the expected observed temper- 
ature T = Ta{r±/rs)~'^^^ is approximately 100 keV due to 
adiabatic expansion. We note that this is a factor of a few 
below the measured peak energy values. However, the as- 
sumption of adiabatic expansion is not necessarily valid since 
energy dissipation occurs below the photosphere to create 
the pairs. As discussed in §3.1 the heating is most probably 
continuous due to oblique shocks. 

Assuming that the photospheric emission that we are 
observing is due to the pair photosphere thus results in a 
value of ro which is similar to the generally assumed value for 
the jet bounding radius. It naturally alleviates the concern 
that the value, found for the baryonic photosphere (ro ~ 10* 
cm), is too large. It can thus be argued that it is indeed 
the pair photosphere we are observing. We note, however, 
that Zhang & Woosley (2004) found typical values close to 
ro ~ 10* cm, since the jet is not well collimated at the center 
thus preventing the acceleration. This value is similar to 
what we deduce for the baryonic photosphere. 



4.2 Epoch 2 

During epoch 2 strong dissipation has to occur at a moderate 
optical depth, r ~ few, in order to broaden the MeV bump. 
From simulating various dissipation scenarios (described in 
§3.3), we conclude, in particular, that the low-energy slope 
is mainly determined by the contribution of synchrotron 
emission (es): For instance, a power-law distribution of 
electrons, produced in the dissipation process, is expected 
to have a peak energy at 7m ~ e^{mp/me) — ISGOte. 
For the typical parameters for GRRB090902B this trans- 
lates into a peak of the synchrotron spectrum at E^"'^'^ — 
186 I/g3^eJ{^j^ee,-iJ?i2^ keV, which is less than the averages 
thermal peak lying at £th<=""''i = 2.82 < kT >~ 840 keV. 
This is still the case if et and are somewhat larger, ap- 
proaching the equipartition values. 

The right-hand panel in Figure [6] shows an example 
spectrum from a simulation of a shock dissipation at the 
same optical depth as before (r — 10), but with an increased 
energy dissipation. In this example, the dissipation fractions 
are given by td = 0.2, te = 0.3, es = 0.3. This spectrum is 
similar in shape to the blue spectrum in Figure[2] illustrating 
that subphotospheric heating indeed can explain the change 
in the spectral shape that is observed in epoch 2. 



4.3 Transition between the two epochs 

Variations at the base of the flow are expected to lead to 
rapid variation of F and rph of the photosphere, down to 
a time-scale of ro^j/c (Rees & Meszaros 2005). This causes 
varying properties of the dissipation, such as its strength 
and where it occurs relative to the photosphere. 

Due to strong dependence on the Lorentz factor of the 
ratio Tsh/fph oc F^ a change in F (apart from variation in 
the other parameters) can easily alter the dissipation pattern 
and thereby the appearance of the photospheric spectrum. 




Ep as re- 
lated to the temperature of the photospliere, this can be explained 
by the expected relation between temperature and spectral width 
in the presence of subphotospheric heating. 



A decrease in F is, for instance, able to cause Tph to be- 
come larger than Vsh, that is, shocks mainly appear below 
the baryonic photosphere - subphotospheric shocks. Such a 
change is hence a plausible explanation the change in the 
spectrum at ~ 12.5 s in GRB090902B. Indeed, in both the 



baryonic (S 4.1.1 1 and pair scenarios (S 4.1.2 1 the bulk of the 



dissipation site then moves from being above rph to being 
below it, mainly due to a drop in F. 

In the photospheric model, the peak of the spectrum 
(which is measured by Ep) should be closely related to the 
temperature of the photosphere. This is, for instance, illus- 
trated by Figure[6] Furthermore, in deriving the Lorentz fac- 
tor, F, from the observables, we note that it most strongly 
depends on the temperature (F oc T^^^, Pe'er et al. 2007): 
The evolution of F and kT are expected to track each other. 
We therefore argue that variations in Ep are closely related 
to corresponding variations in F (apart from the Doppler 
boost). In Figure |4] we show that there is a drop in averaged 
-Bp between epoch 1 and epoch 2. This thus indicates that a 
drop occurred in the averaged temperature at this time and 
thereby also a drop in F. 

The change in F leads to a change in dissipation pattern, 
according to the discussion above, and thereby the spectral 
shape. Indeed, Figure [7] shows a correlation between Ep and 
a, reinforcing this interpretation; a decrease in Ep corre- 
sponds to a decrease in F, which leads to more subphoto- 
spheric heating thereby broadening the spectrum (a decrease 
of a) . This is most clearly illustrated by the local increase in 
Ep observed around 15 s (see Fig.|4|, which corresponds to a 
local increase in a, indicating less subphotospheric heating. 



5 DISCUSSION 

The remarkable property of GRB 090902B is that it is dom- 
inated, at early times, by a very peaked spectral component 
in addition to a power law component. This peaked com- 
ponent, which has a spectral shape very close to a Planck 
function, can only be explained by emission from the photo- 
sphere. As we have seen in the analysis above, at late times 
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Figure 6. Generic model spectra from subphotospheric shock heating (solid line), illustrating the broadening effect. Left panel, A low 
level of dissipation only slightly distorts the MeV peak from begin a Planck function. This spectrum can explain the observed shape 
at 6.5 s (see Fig. [2|. Right Panel. A higher level of dissipation broadens the photospheric peak, leading to a Band-like spectrum. This 
spectrum can explain the spectrum at 22 s (see Fig.[2]|. The dashed, red line shows the shape of a Planck spectrum at the temperature 
that corresponds to the peak of the spectrum. 



the spectrum of this peaked component is broadened into a 
spectrum that can be described by a Band function with an 
a slope of approximately -0.6 and a P slope of approximately 
-2.5. We note that this is close to the spectral shape usually 
attributed to GRBs (Kaneko et al. 2006). Since there is a 
continuous change in shape, even though it occurs over a 
short period compared to the burst duration, it has to be 
concluded that the photosphere emission continues during 
the second half of the burst. The consequence of this con- 
clusion is that the photosphere emission spectrum can have 
a large variety of shapes, which reflect the dissipation pro- 
cesses in the flow. We have argued that this burst provides 
observational evidence for subphotospheric heating. 

This mechanism provides a natural explanation to the 
observed variety of spectral shapes in GRBs. In addition, the 
more typical spectral evolution, in which the spectra grad- 
ually become softer, can be explained by a gradual change 
in the dissipation pattern in the flow. 

The shape of the photospheric peak, given, for instance, 
by Q and /3 of a Band function fit, can now be translated 
into physical properties of the dissipation, quantified by the 
parameters such as e^, e_g, and r and dissipation rates. De- 
pending on the details of the model we can therefore diag- 
nose the outflow and its temporal evolution in individual 
bursts. 



5.1 Photospheric emission 

Among the main contenders for explaining the hard a val- 
ues observed are (i) small-pitch-angle synchrotron emission 
(Epstein 1973; Epstein & Petrosian 1973) or similarly jitter 
radiation (Medvedev 2000, 2006), and (ii) inverse Comp- 
ton emission seeded by self-absorbed synchrotron emission 
(Painatescu & Meszaros 2000) or by soft photons with a nar- 
row energy distribution, i.e., a quasi-monoenergetic distribu- 
tion (Stern & Poutanen 2004). However, these non-thermal 
emission models typically lead to very broad spectral peaks. 



and cannot produce spectra that are as narrow as observed 
(below an order of magnitude) and that are as hard (a ^ 0, 
see [2.2\. For instance, inverse Compton emission leads to 



broad spectral peaks of typically 2 orders of magnitude or 
above (see, e.g. Stern & Poutanen 2004, Baring & Braby 
2004). 

Moreover, studies of the acceleration processes in rel- 
ativistic collisionless shocks indicate that a strong ther- 
mal component in the electron spectrum is formed (Sironi 
& Spitkovski 2010). Even though the full emitting region 
can still not be fully simulated, these results question the 
classical assumption of a strong power-law component of 
the shocked electron population. Sironi & Spitkovski (2010) 
show through particle-in-cell simulations of shocks in un- 
magnetised pair plasmas that the width of the emitted syn- 
chrotron spectrum, in general, lie in the range of 2 - 4 orders 
of magnitude. 

We further note that Sironi & Spitkovski (2010) argue 
that if the electrons are accelerated in relativistic unmagne- 
tised shocks then the emission is in the classical synchrotron 
regime rather than in the jitter radiation regime. In addition, 
small-pitch-angle synchrotron emission predicts a negative 
correlation between Ep and a, in contrast to the observed 
(mainly) positive correlations (Lloyd-Ronning & Petrosian 
2002). These two facts pose further challenges for such an 
explanation for the hard spectra observed. 

The narrow and hard spectra we observe in 
GRB090902B are thus inconsistent with the non-thermal 
emission models. 

The photospheric model can easily overcome many of 
the challenges of the standard, internal-shock synchrotron 
model (see, e.g., discussion in Ryde & Pe'er 2009). Most im- 
portantly, the (reprocessed) Planck function naturally pro- 
vides very hard spectral slopes (a ^ 1). Moreover, a power- 
law distribution of electrons is not required to be produced 
by the acceleration processes; the power-law slope in the 
BB+pl model has a preferred value of a ~ —1.5 (Battelino et 
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al. 2002, Ryde & Pe'er 2009). This is naturally expected due 
to the cooling of the electrons, which produces a power-law 
distribution below the characteristic synchrotron frequency 
(the case is similar for SSC). 

An inevitable and characteristic signature of emis- 
sion from the photosphere is a cut-off at energies 
O.STMeV ^ 100 MeV due to pair production. The obser- 
vation of an extension of a Band function to energies above 
several GeV in bursts (e.g. GRB080916c, Abdo et al. 2009) 
may thus pose a challenge for (reprocessed) photospheric 
emission to explain the spectra. In order to overcome this 
difficulty, the existence of very high-energy photons in the 
observed spectra therefore requires an additional emission 
site, which is capable of producing optically-thin emission. 
We note that such a scenario (spectra indicating two emis- 
sion components) is clearly observed in several bursts ob- 
served by Fermi-LAT. Moreover, in one of the analysed time- 
bins in GRB080916c there is indeed an indication of an extra 
power-law component, even though the overall conclusion of 
the analysis is that a single component sufficiently fits that 
data (Abdo et al. 2009). This suggests the possibility that an 
extra component could indeed be present in the data. How- 
ever, the combined spectrum can still be satisfactorily fitted 
by an extension of a single Band function, found at lower en- 
ergies. Ghisellini et al. (2010) indeed find that the power-law 
slope fitted to the spectrum at energies ^ 40 MeV (i.e., the 
Band /3) and the power-law slope fitted to the spectrum at 
energies at ^ 100 MeV, are significantly different. This can 
be an additional indication that additional emission compo- 
nents are required at the highest energies, thereby alleviat- 
ing this problem for the photospheric model. 



5.2 Correlation between Ep and a 

In j ]4.3| we discuss the observed correlation between Ep and 
a in terms of subphotospheric heating. A positive correla- 
tion between time-resolved Ep and a values for individual 
bursts was early identified by, for instance. Ford et al. (1994), 
Crider et al. (1997), and Lloyd-Ronning & Petrosian (2002). 
More recently, Kaneko et al. (2006) found that this is the 
strongest correlation among GRB parameters. In a sample 
of 196 bursts they found a strong correlation for 26 %; in 
most cases the correlation is stronger than the one in Figure 
[T] The actual fraction of bursts exhibiting a strong correla- 
tion can be even higher, since measurement uncertainties in 
many cases may have masked the correlation. 

The fact that not all bursts exhibit a strong correlation 
between Ep and a might also be an indication of that the 
observed spectral peaks are not directly related to kT. In 
some bursts a weak contribution of a Planck function on 
top of a dominant Band spectrum can be identified (Ryde 
& Pe'er 2009, Guiriec et al. 2011). The power peaks are, 
in these cases, given by the non-thermal Band function and 
should thus not be identified by kT. As discussed in §3.2| 
this could be the case if the energy density of the electrons is 
much larger than the energy density of the thermal photons. 
Then the thermal component is expected to be relatively 
weak. 



5.3 Change in spectral shape from "thermal" to 
" non-thermal" 

Other scenarios than the one described in §4.3| can be envi- 
sioned to describe a change from quasi-Planck spectrum to 
Band-type spectrum. For instance, Beloborodov (2010) cal- 
culated the full radiative transfer of the relativistic jet and 
showed that a perfect Planck function is obtained if and only 
if the flow is dominated by radiation, that is, the radiation 
energy density is much larger than the rest mass energy den- 
sity of the plasma in the comoving frame of the flow. This is 
independent of subphotospheric dissipation since the photon 
energy totally dominates. However, there will inevitably be 
broadening due to geometrical effects (Pe'er & Ryde 2011), 
as discussed in §3.2[ The spectral change observed between 
epoch 1 and 2 can then be due to a change from a radiation 
dominated phase into a baryon dominated phase, in which 
dissipation causes the strong deviation from the Planck func- 
tion. An interesting consequence of such a scenario would 
be vanishing of the strong polarisation expected during the 
photon dominated phase. 

Yet another possibility is if the variability time-scale 
and amplitude suddenly change. During epoch 2 the varia- 
tions of the temperature can become large and be on a time 
scale smaller than the integration time scale. We would then 
measure spectra that are broader than a Planck spectrum 
due to significant variations of the spectral peak during the 
integration time. In such a case the asymptotic slopes of the 
spectrum should still be those of a blackbody. 

Finally, one may still envision that the emission in 
the MeV peak during epoch 1 is due to the photosphere 
while during epoch 2 the emission is from a different, non- 
thermal, radiation process, for instance, from optically-thin 
synchrotron emission. However, as discussed above the ob- 
served spectra during epoch 2 are still significantly inconsis- 
tent with what is expected from the fast cooling electrons 
(a — —1.5 and a much broader spectral width). Moreover, 
to get the synchrotron peak-energy to lie in a similar en- 
ergy range as the thermal peak during epoch 1 requires an 
unreasonable coincidence. This is also not supported by the 
behavior of the NT part of the spectrum. 



5.4 Comparison with Pre-Fermi Spectral Analysis 

Prior to the launch of Fermi we mainly had to rely on spec- 
tra in a narrow energy range, for instance that of CGRO- 
BATSE (~ 25-1900 keV). The limited spectral widths made 
it difficult to unambiguously deduce the spectral behaviour 
of GRBs. Nevertheless, the importance of the photospheric 
emission was already alluded to (see, e.g., the review by 
Ryde 2008). Using the broader energy range provided by 
Fermi this unambiguity can now be revealed. We can sum- 
marise the pre-Fermi results with the following three be- 
haviours. 

(i) Many bursts provided us with a clear indication of 
the photosphere. These include the spectra that are well de- 
scribed by a single Planck function (Ryde 2004) and the 
spectra for which a fit of the BB-|-pl model gives a statisti- 
cally significant improvement over a fit with a Band function 
model (Ryde & Pe'er 2009). Such cases should be similar to 
GRB090902B observed by Fermi; the vF^, peak is due to the 
Planck function and there are two distinct components. 
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(ii) Other bursts indicated that the power-peak lay be- 
yond the observed energy range, even though a subdomi- 
nant, thermal peak was identifiable: The thermal component 
forms a shoulder on the low-energy side of the power-peak. 
The power peak of the spectrum is, in this interpretation, 
not directly due to the photospheric component, but due to a 
non-thermal emission. Example of such spectra are given by 
PHEBUS Granat observations of GRB 900520a (Barat et al. 
1998), CG7?0 BATSE/EGRET observations of GRB981021; 
BATSE trigger 7071 (Gonzalez et al. 2009; Ryde & Pe'er 
2009), and Fermi observations of GRB100724B (Guiriec et 
al. 2011), as well as Fermi cases like GRB 080916c (Abdo et 
al. 2009). 

(iii) As argued in the study above, the thermal peak 
can be broadened due to subphotospheric heating, which 
creates a spectrum which is Band-like and significantly dif- 
ferent from a Planck function. This mechanism is thus able 
to explain many Band-like spectra as well as the typical 
hard-to-soft spectral evolution. Apart from these facts, it 
also opens up the possibility that some spectra that were 
fitted in pre-Fermi analyses by a BB-|-pl model, in a narrow 
energy range, may have been misinterpreted. This could be 
the case particularly for bursts whose BB-|-pl spectral fits are 
not an improvement over a Band-only model fit. In these 
cases, the Planck function component of the fits still cor- 
rectly captures the photospheric peak. However, the power- 
law component instead captures the broadening of the peak 
(which makes it non-Planckian), instead of representing a 
real secondary component. If this is the case, such a BB-|-pl 
model would, of course, not be able to describe the broad 
band spectrum of these bursts. Such spectra can explain 
the results in Ghirlanda et al. (2007) and Bellm (2009), who 
find that in some bursts a simple extension of the BB-|-pl 
model does not seem to fit the data at hand outside the 
BATSE energy range. We note that this is apart from the 
fact that the non-thermal component is expected to be a 
more complicated function than a single power-law over a 
broader energy range. 



6 CONCLUSION 

The unambiguous signature of emission from the photo- 
sphere is a Planck function. However, subphotospheric dis- 
sipation can easily distort the photospheric emission into a 
broader spectrum, resembling a Band function. 

The burst of 090902B made it possible for us to draw 
two important conclusions. First, the study of the behaviour 
of the MeV peak allows us to observationally claim that the 
origin of a Band function can be the same as the origin 
of a Planck spectrum; a Band spectrum can be interpreted 
as photospheric emission. Second, in order to explain the 
broader spectral shape of the photospheric emission, sub- 
photospheric dissipation is needed. The existence of such 
dissipation is thus verified. 

We show that a varying Lorentz factor of the outflow 
leads to a varying strength of subphotospheric dissipation. 
This, in its turn, leads naturally to a correlation between 
the broadening of the spectrum and its peak energy. A cor- 
relation between a and Ep, which is a consequence of this, 
is indeed commonly observed in GRBs. 

Photospheric emission could thus be a ubiquitous sig- 



nature of the prompt emission spectrum and play an impor- 
tant role in creating the diverse spectral shapes and spectral 
evolutions that are observed. The photospheric component 
can be dominant (e.g. GRB090902B) or subdominant (e.g. 
GRB100724B). It can be a Planck function throughout the 
burst (Ryde 2004) or it can broaden with time (as argued 
in this paper). An important consequence of this is that 
broad-band fitting of GRB spectra with only a Band func- 
tion might be misleading. 
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